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1 Introduction 
1.1 

This paper is about a certain interaction between probability and geometry. 
The random objects involved will be random stepped surfaces spanning a 
given boundary in R 3 . Equivalently, one can talk about random rhombi 
tiling of a planar domain or random dimer coverings of certain subgraphs of 
the hexagonal graph. Probabilistic questions about these random surfaces 
will be answered in terms of a nonrandom algebraic object, geometrically a 
curve in a noncommutative plane. 

Underlying this connection is Kasteleyn's theory of planar dimers which 
computes all probabilities in terms of the Green's function of a certain finite- 
difference operator K. It will be clear from our construction that the connec- 
tion between finite-difference operators and noncommutative geometry may 
be easily extended far beyond what we do in this paper. Our goal here, how- 
ever, is to explain a certain phenomenon in the least possible generality and 
to stay as close as possible to certain specific applications. These applica- 
tions, as well as some other directions that look promising will be discussed 
below. 



1.2 

Let f2 be a simply-connected planar domain that can be tiled by rhombi as 
shown in Figure [U Well-known bijections illustrated in Figured] identify such 
tilings with dimer coverings of subgraph fig = fi PI of the hexagonal graph 
r$ and also with stepped surfaces spanning given boundary. An introduction 
to dimers and stepped surfaces may be found in [8|. 
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Figure 1: Stepped surfaces are the same as tilings of a planar domain Q by 
rhombi. Each tile is a union of a black and white triangle. The adjacency graph 
of the triangles is (a piece of) the 6-gonal graph. 

Stepped surfaces arise in mathematical physics in a variety of contexts: 
from simply-minded, but realistic models of interfaces (e.g. crystalline sur- 
faces) to the sophisticated setting of super-symmetric gauge and string the- 
ories (see e.g. [TU [16] for an introduction). 

In all applications, it is natural to weight the probability of a stepped 
surface S by the volume V(S) enclosed by it, i.e. to set 

Prob(S) oc q vis) , (1) 

where q > is a parameter. Note that for two surfaces Si and S 2 spanning 
the same boundary the difference V(Si) — V(S , 2) is well-defined, which is all 
that matters in (CQ). In the crystal surface context, logg is the energy price 
for removing an atom0 

1.3 

We will be particularly interested in the case when Q grows to infinity, while 
keeping its shape, that is, the number and orientation of its boundary seg- 
ments. We will refer to such domains as polygonal. The study of a random 
stepped surface a large polygonal boundary (equivalently, random dimer cov- 
ering, or random tiling) leads to interesting probabilistic questions. The 
nontriviality of these questions may be appreciated by looking at Figure El 

1 Note that there is a well-known ambiguity in reconstructing 3-dimensional surfaces 
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Figure 2: A simulation of the limit shape formation. The curve separating the 
ordered regions (facets) from the disorder is called the frozen boundary. Here, it 
is a cardioid, i.e. the dual of a rational cubic. 

Apparent in Figure [2] is a formation of a certain nonrandom limit shape, 
in other words, as the mesh size goes to zero so does the scale of randomness. 
This is a form of the law of large numbers. The existence of a limit shape was 
proven for stepped surfaces (with arbitrary boundary conditions) by H. Cohn, 
R. Kenyon, and J. Propp in [3]. 

In [9], the limit shape for polygonal boundaries was linked to a certain 
plane algebraic curve Q. In particular, the frozen boundary, which is the 
curve separating order from disorder in Figure [21 is the planar dual of Q in 
exponential coordinates. 

1.4 

The main object of this paper may be characterized, informally, as a quanti- 
zation of the limit shape, or of the curve Q to be more specific. This quanti- 
zation exists for finite Q, that is, before any limits are taken. In particular, 
it captures not just the limit shape but also the fluctuations of our random 

from tiling, namely, one can switch the roles of convex and concave corners. A rotation 
by 7r/3 interchanges the choices. Our conventions are fixed by (f2]). 
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surfaces. Viewed like this, it should not be surprising to see noncommutative 
curves appear. Note, however, that technically the noncommutativity will 
be linked to the parameter log q and not to the size of fluctuations (as could 
be expected from the uncertainty principle). 

1.5 

In principle, classical Kasteleyn's theory [6] answers all possible questions 
about random stepped surfaces in terms of the Green's function, i.e. the in- 
verse of a certain difference operator. This Kasteleyn operator K is a weighted 
adjacency matrix of the graph fi 6 . 

Note that r 6 , and hence, Q e is bipartite, that is, its vertices may be 
colored in two colors (black and white, traditionally) so that only vertices 
of opposite color are joined by an edge. This is reflected in Figure [31 We 
will index the rows (resp. columns) of the adjacency matrix by white (resp. 
black) vertices. The nonzero matrix elements Kij should satisfy 

gK 21 K 43 K 65 = K 23 K 45 K 61 (2) 

for each face F, where Vi, . . . ,Vq are the six vertices going around F as in 
Figure [31 This fixes K uniquely up to a certain gauge transformation, namely 
left and right multiplication by a diagonal matrix. 




Figure 3: Each hexagon face carries logq units of magnetic flux. 

The goal of this paper may be informally described as looking for some 
hidden structures in the inverse matrix K -1 . Certain structures in K _1 are 
plain to see: by definition, the entries of K _1 satisfy a finite-difference equa- 
tion in each index, namely K K _1 = K _1 K = 1. 
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Our main claim is that for polygonal domains Q the entries of K _1 satisfy 
additional finite-difference equations. The degree of these additional equa- 
tions is determined by the shape of Q, that is, by the number of boundary 
segments, and not by the size of Q. This is crucial from the probabilistic 
viewpoint. 

1.6 

The noncommutative geometry of the title provides a natural language to 
state and study these additional equations. 

The origin of the noncommutativity may be traced to (J2J). For q = 
1, the adjacency matrix of T 6 is an obvious solution and this solution is 
translation- invariant, i.e. commutes the the subgroup Z 2 C Aut^) acting 
by (bipartition-preserving) translations. 

For 5 / 1, the translation- invariant equation (j2J) has no translation- 
invariant solutions, which means that the Kasteleyn operator K now com- 
mutes with magnetic translations, i.e. translations followed by a gauge trans- 
formation. In turn, magnetic translations commute only up to a factor 
(whose logarithm is proportional to the area of the parallelogram spanned 
by the translation vectors). They form, in other words, an algebra known as 
the quantum 2-torus. 

The considerations so far apply only to the whole 6-gonal graph T 6 , that 
is, in the absence of any boundaries. Somewhat remarkably, however, a cer- 
tain framework may be established in which the Kasteleyn operator and the 
commuting magnetic translations act in a way compatible with polygonal 
boundaries. This involves compactifying the quantum 2-torus to a noncom- 
mutative plane, meaning that one introduces a certain graded algebra A, a 
deformation of the ring of polynomials in xi,X2,Xs, such that the quantum 
2-torus is the degree part of A [( y xiX2Xs)^ 1 ]. 

1.7 

For any fixed white vertex w G fi 6 , the action of magnetic translation on 
K" x (- ,w), that is, on the corresponding column of K , yields a graded A- 
module Q w . The additional equations satisfied by K _1 will be reflected in the 
fact that Q w is a torsion module. 

For different w, the modules Q w share the same fundamental features. In 
fact, there is a canonical submodule Q in all of them that depends on Q only 
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and captures the essential information. 
1.8 

The construction of the module Q and the study of its basic properties will 
occupy the bulk of the present paper. While the definition of Q involves 
nothing beyond elementary combinatorics and linear algebra, we will find 
the resulting object has a certain depth and complexity. 

The degrees of its generators and relations (and hence the degrees of the 
additional equations satisfied by K _1 ) are determined by the combinatorics 
of the domain Q only, see Theorem [TJ On the other hand, the explicit form 
of these relations depends of q and the geometry of Q in a rather intricate 
fashion. 

1.9 

This is where a geometric way of thinking about such modules, pioneered 
by M. Artin and his collaborators, becomes essential (see e.g. [21] for an 
introduction). While its a matter of definitions to associate to Q a sheaf 
on the noncommutative plane, the geometric intuition thus gained is very 
valuable. 

In the first place, this is what allows us to view Q as a quantization of the 
limit shape Q. In fact, quantization requires additional degrees of freedom, 
parametrized by line bundles (and more general rank 1 sheaves) on Q. These 
may be compared to complex phases in quantum mechanics. 

Certain features of Q, such as its points of intersection with the coordinate 
axes of P 2 have a direct quantum analog satisfied by Q, see Section |U Others 
(such as the geometric genus) are harder to generalize, see [TO] . 

1.10 

Because Q is a purely combinatorial object, we can modify it by simply 
moving boundary segments in and out. We prove in Section that these 
transformations act on Q by what may be called a noncommutative shift on 
the Jacobian. In particular, this describes what happens to Q under rescaling 

of a 
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2 The quantum limit shape 
2.1 

Let A be the algebra generated by Xi,x 2 ,x 3 subject to the relations 

XjXi Qij XiXj , (3) 

where q ti = 1 and qv,- = qj t . This is a basic example of a noncommutative 
projective plane P 2 , see e.g. [21] for an introduction. Adding the inverses x~ x 
of the generators, we obtain a larger algebra T known as a noncommutative 
3-torus. 

Both A and T are graded by the total degree in all three generators. 
Let Td, where d G Z is arbitrary, be one of the graded components. The 
monomials 

™a _ ^ a 2 a 3 _ -r 

obviously correspond to lattice points a G Z 3 lying in the plane a\ + a 2 + 
<23 = d. Their projection in the (1, 1, 1) direction, that is, their image in 
R 3 /R ■ (1, 1, 1) = M 2 may be identified with the black vertices of T 6 . In the 
same fashion, the monomials in T^+i are put into bijection with the white 
vertices of r 6 . This is illustrated for d = 4 in Figure |U 
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Figure 4: Black and white vertices may be identified with points of Z 3 lying on 
two parallel hyperplanes and thus with monomials of two successive degrees. Here 
we only plot points in R^q which correspond to polynomials in the Xj's. 

2.2 

Consider the operator 

K : T d — > Jd+i 

given by right multiplication by x\ + X2 + x-s G Wi, that is, 

K/ = / • (sj + ar 2 + ar 3 ) • 

One easily checks from the commutation relations (j3J) that, indeed, in the 
basis of monomials, this is a g-weighted Kasteleyn operator for T 6 as above, 
with 

q = qi2 <?23 931 • (4) 

Note that while there is no canonical ordering of the variables and, hence, 
no canonical normalization of a monomial, the lines spanned by monomials 
in T are well defined. This is all we need since we only care about K modulo 
gauge transformations. 

2.3 The module M 
2.3.1 

Monomials are normal in T, that is, Ax c = x c A is naturally a A-bimodule. 
Assuming degx c < 1, monomials in (Ax c )i form a triangle with vertices 

™Cl C2 -C1-C2 + 1 Ci -C1-C3 + I C3 „— C2— C3+1 C2_.C3 g- T 
1 2 3 ? 2 2 3 ' 1 2 3 1 * 
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We denote this equilateral triangle by A(c) and call it the support of (Ax c )i . 
2.3.2 

We think of black and white vertices in f2 6 as monomials in T and Ti, 
respectively. Write Q as a set-theoretic combination of triangles 

Q = |jA(a,;)\|jA(b J ), a^b.-eZ 3 , (5) 

* i 

as in Figure EJ The inclusion of triangles induces the inclusion of graded 
A-bimodules 

0Ax b 'C0Ax a '. (6) 

We denote by M the quotient A-bimodule in (jf)]). Define an endomorphism 
K of M as the right multiplication by x\ + x<i + x% G Ai . This is a map of 
left A-modules. By construction, the operator 

K : M -> Mi 

is the Kasteleyn operator for Qq. 



Figure 5: The domain f2 may be represented as a difference of unions of triangles. 
Here, it is the union of 2 blue triangles minus the union of 4 green ones. 
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2.3.3 



The idea of stable range will play an important part in this paper. By 
definition, a number i is in the stable range for M if the domain 

fl(i) = supp M i+ i 

has the same combinatorics as Q = Q(0). In particular, the length of the 
shortest white boundary of Q, i.e. a boundary formed by white triangles, 
gives an upper bound on the stable range. 

More generally, refining ([5]), M may be given a resolution 

F 3 F 2 -> F 1 -> F -> M -> (7) 

by modules of the form ® A , where the maps are the natural inclusions. 
Here Fq = 0Ai a ' corresponds to the triangles A(aj) in ((Sj). The intersec- 
tions among Ai a ' together with A(bj) contribute to relations Fx, and so on. 
The module F doesn't have generators of positive degree, but the other Fj's 
do. The stable range is until the first such generator appears. 

Note that the stable range scales linearly with Q, meaning that it scales 
like the inverse mesh size in the probabilistic setup. Throughout the paper, 
we will only be interested in what happens in the stable range. 



2.3.4 

Our next goal is to compute the Hilbert function of in the stable range. 
For that, we will make a genericity assumption that all 3 boundary slopes 
appear in the boundary dVL of Q in a cyclic order. We will denote by deg Q the 
number of times dfl cycles through the three slopes. For example, deg Q = 3 
in Figure [61 

Lemma 1. For al > in the stable range, we have 

dimM d = dimM - degf2 ^ 2 - - d ind K . 

Here, of course, the index of K equals zero, but it will not vanish in the 
generalizations considered below. 

Proof. Consider the map 

M d 3f^ x 3 f E M d+1 . (8) 
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The kernel and cokernel of this map are formed by functions supported on 
horizontal strips of the form shown in Figure O More precisely, white hori- 
zontal boundaries correspond to the cokernel, the black ones — to the kernel. 




Figure 6: A white (left) and black (right) boundary strips. Note that they always 
have the same shape provided the boundary of ft, cycles through the 3 directions. 

Note that white boundaries of supp shrink with i, while the black ones 
expand. Thus each of the degO horizontal boundaries contributes —1 to the 
second derivative of dimMj with respect to i. We conclude 

d 2 

dimM d = - degQ— + ... , . 

Here dots stand for a polynomial in d of degree < 1, which is uniquely fixed 
by its evaluation at d — 0, 1. □ 

2.3.5 

Lemma 2. The operator K : Md — > Md+i is surjective for d > in the stable 
range and q > or q generic. 

Proof. For d = 0, by Kasteleyn's theorem, det K determinant gives the q- 
weighted count of stepped surfaces, hence nonzero for q > 0. For d > 0, 
we use the commuting map ([HD and induction. By induction, it suffices to 
check the surjectivity of K on a white boundary strip, as in Figure El which 
is immediate. □ 

2.4 The module Q and the inverse Kasteleyn matrix 
2.4.1 

We denote by Q the kernel of K acting on M. This is a graded left A-module. 
To see why this may be a useful definition, let us generalize the construction 
slightly. 
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2.4.2 



Let w be a white vertex of corresponding to a monomial x w . Let f2 w be 
obtained from Q by removing the corresponding white triangle and set 

M w = M/Ax w . 

The conclusions of Lemmas [1] and [2] continue to hold for M w , with the obvious 
modification that 

degfT = degfi + 1, indK£ = l. 

In contrast to Q, we have Q™ ^ 0. Indeed, by construction, Q™ is spanned 
by the corresponding column K~ x x w of the inverse Kasteleyn matrix. 

Via the left A- module structure on Q w , the algebra A acts on the columns 
of the inverse Kasteleyn matrix by difference operators. To see that a nonzero 
difference operator from A must annihilate Q™ , it will suffice to compute the 
Hilbert function of Q w . 

2.4.3 

From Lemmas [1] and [2] we have, in the stable range, 

dim Q d = d deg Vl + ind K = d deg f2 (9) 

and, similarly, 

dimQ£ = (degfi + l)d + l. (10) 
Since this dimensions grow only linearly in d, for any g e Q d the map 

Ai 9 / i-> f-ge Q d+i 

must have a kernel as soon as i is large enough. These are the sought differ- 
ence equations satisfied by K -1 . We will, obviously, have do more work to 
say something more specific about them. 

2.4.4 

Note that we have an exact sequence 

O^Q^CT^L^O, (11) 
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where the third term L satisfies dimL d = d + 1. In fact, L is a line module, 
i.e. the module of the form 

L = A/M, ZeAi. 

Lemma 3. In the stable range, i.e. for w sufficiently far from the boundary 
of Q, L is a line module with I = x w {x\ + x 2 + £3) £~ w . 

Proof. Let g = K -1 :r w be the generator of Qq and suppose fg e Q for some 
/ e A. This means that may be extended to all of as a solution of 
Kasteleyn's equation. In other words, there exists a polynomial g' G Ad eg /-i 
such that 

= Hfg + g'x w ) = fx- + g'x w ( Xl + x 2 + x 3 ) , 
proving the assertion. □ 

From this perspective, there isn't much difference between Q and Q w . 

Somewhat poetically, we will call Q the quantum limit shape. Mathemat- 
ical reasons for this name will be discussed below. 



3 The structure of Q 
3.1 

The goal of this Section is to prove the following 

Theorem 1. In the stable range and for generic q, Q is generated by degf2 
generators of degree 1 subject to deg Q linear relations. In other words, the 
minimal graded free resolution of Q has the form 

-> A(-2) dcgn -> A(-l) dcgn -> Q -> . (12) 

Similarly, 

-> A(-2) dcgn -> A © A(-l) dcgn - 1 -> Q w -> . (13) 
Here we denote, as customary, A(i)^ = Aj + d. 



13 



3.2 



In the commutative case, resolutions of the form ffl2|) are well known in 
algebraic geometry, see, in particular, pQ. The corresponding sheaves on P 2 
are of the form i*£, where 

is an inclusion of a curve of degree D = degf2 and £ is a line bundle £ (or 
a more general torsion-free sheaf in case Q is singular) of degree g — 1 . Here 

0=(Z>-l)(£>-2)/2 

is the arithmetic genus of Q- Concretely, 

Q = det R , 

where R is the matrix of linear forms that gives the map 

R : A(-2) D -> A(-l) D 

in (112j) . The condition in ( fl2l) on £ to have no sections means £ G Jac s _i(Q)\ 
9, where 9 C Jac 9 _i(Q) is the theta divisor. 

The meaning of ( Tl~3l) is parallel, with the difference that D becomes 
degf2 + 1 and deg£ becomes g. In the commutative case, (ITTj) implies the 
support of Q w has the line I = as a component. 

3.3 

Cne of the eventual goals of the present project is to understand the behavior 
of Q as the mesh size goes to while logg — > at a comparable rate. The 
moduli of sheaves of the form fTl~2l) have a natural compactification, which 
guarantees that any sequence has a subsequence converging to an actual sheaf 
on the usual commutative plane P 2 . 

It will be shown in [10] that this limit is supported on the curve Q cor- 
responding to the limit shape. This is the reason for calling Q the quantum 
limit shape. 

3.4 

For any graded A-module Q, the degrees of its generator, relations, etc., may 
be read off the dimensions of the graded components of the vector spaces 

Tor J Q = Tor l (A/A >0 ,Q), 
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where A >0 is the ideal generated by Ai. In particular, the existence of a free 
resolution of length 2 is equivalent to the following 

Lemma 4. In the stable range of degrees, Torj Q = for % > 1. 

Proof. Since Tor>3 vanish identically for the algebra A, we need to show the 
vanishing for % = 2, 3. Since the resolution © is defined combinatorially in 
terms of Q, we conclude Torj M = in the stable range. This is, really, the 
definition of the stable range. Since Q C M and Tor 4 = identically, we 
conclude 

Tor 3 Q = . 

By construction, there is an exact sequence 

and since ImK C M, we have Tor3 Im K = 0. Now from the long exact 
sequence for Torj, it follows that Tor 2 Q = 0. □ 

3.5 

Lemma 5. For generic q, Q is generated by Qi. 

Proof. It suffices to consider the commutative 1 . Then, on the one 

hand, Q is annihilated by x\ + x 2 + x$, while on the other Torj Q = 0, i > 1. 
This means Q corresponds to a vector bundle on the line x\ + x 2 + £3 = 0. 
From its Hilbert polynomial, we see that it must be 0(— l) degn , whence the 
conclusion. □ 

3.6 

Now it is easy to complete the proof of (1121) . We have degf2 = dimQi 
generators in degree 1 and from ([9]) we see that they must satisfy degf2 
linear relations. There are no other generators by Lemma [5] and no other 
relations by (jUJ). 

The proof of (fT3]) goes along the same lines. Lemma H] still holds even 
though (ToriM w )i ^ 0. The analog of Lemma is that Q w is generated 
in degrees and 1, because in the commutative case Q corresponds to the 
bundle © 0(-l) dc s Q on the line x 1 + x 2 + x 3 = 0. 
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It remains to explain why the commutative resolution 

-> A(-l) © A(-2) dcgf7 -> A © A(-l) dogn -> Q w -> , g# = 1 , 

jumps to f)13p for generic g. In other words, we need to check that the 
generator of Qq no longer satisfies a linear relation for generic q. This is an 
easy consequence of the results of the next section. 

Namely, a linear polynomial meets X1X2X3 in 3 points, that is, it annihi- 
lates 3 generators of point modules over A. As we will see, the annihilator of 
Qo meets each coordinate axis in deg 0, + 1 points. 

4 Boundary points 
4.1 

Recall from [9] that the curve Q describing the limit shape is determined 
as the unique rational curve of degree degf2 for which the dual curve Q v is 
inscribed in Q. This means, in particular, that Q meets each coordinate line 
of P 2 in deg Q specified points. 

In this section, we will see that the quantum limit shape Q satisfies the 
exact noncommutative analog of this incidence. 

4.2 

We define 

d 3 Q = Q/x 3 Q . 

We will see that d 3 Q is a direct sum of degf2 point modules labeled by the 
horizontal boundaries of Q as in Figure El 

By definition, the Hilbert polynomial of a point module is equal to the 
constant 1. Up to modules of finite length, point modules are parametrized 
by the toric divisor X1X2X3 = of P 2 . The correspondence is a follows 

(di : a 2 : 0) <-> A/ (x 3 , a 2 Xi — aiX 2 ) ■ 

The ratios a 2 ja\ for the summands of c^Q will be determined by the vertical 
coordinate of the horizontal boundaries of fl. 
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4.3 



The Hilbert function evaluation 

dim(d 3 Q) d = degQ, d > , 
is a consequence of the following 

Lemma 6. Left multiplication by x 3 has no kernel acting on Q. 

Proof. A polynomial in the kernel of left multiplication by X\ has a support 
in a strip of width 1 along the black boundaries, as in Figure El right. It is 
impossible for such function to be annihilated by the Kasteleyn operator. □ 

4.4 White boundaries 

Let Q' be a nontileable domain obtained by moving one of the white horizon- 
tal boundaries of Q one step in. Denote by M' the corresponding monomial 
module. Let K' : M' — ► M' be right multiplication by x\ + x 2 + x 3 and let 
Q' = Q n W be its kernel. Since x 3 Q C M', d 3 Q surjects onto Q/Q'. 

By cutting white boundary strips off fi'(i), i > 1, we see as in the proof 
of Lemma [2] that K' surjects onto M[ for z > 1. Since ind Kq = —1, it follows 
from Lemma [T] that Q/Q' is a point module. 

Lemma 7. Let a be the vertical coordinate of the strip Q \ Q' , that is, let 
deg X3 f = a for all f e M/M'. Then 

Q/ Q' ^ A/ (x 3 , Xi + q a x 2 ) (1), 

starting in degree 1. 
Proof. The series 

{x\ ~t~ X2) — X~[ — X^ X^X-^ ~t~ X-y X 2 X^ X 2 X-^ — ... 

gives 1 after left or right multiplication by (x\ + x 2 ). Interchanging the roles 
of X\ and x 2 and taking the difference 

5 = (xi + x-iY 1 - (xi + x 2 ) _1 

we get an analog of the usual 5-function, which is annihilated by both left 
and right multiplication by {x\ + x 2 ). 



17 



An element of (Q/Q') d is a truncation of the series x^x 1 a+ +1 5 on both 
sides. Left multiplication by 

x x + q a qtt^2 

annihilates it, whence the conclusion. □ 
4.5 Black boundaries 

Now let Q f be a nontileable domain obtained by moving one of the black 
horizontal boundaries of Q one step out. We have a map M' — > M corre- 
sponding to the restriction of functions and from Lemma [6] we conclude that 
it yields an injection Q' — > Q. Further, x 3 Q is in the image of Q', hence Q/Q' 
is again a point module onto which c^Q surjects. 

Lemma 8. Let a be the vertical coordinate of the strip Q' \ Q, that is, let 
deg X3 f = a for all f in the kernel of the restriction map M' — > M. Then 

Q/Q' = A/ (x 3 ,x 1 + q a x 2 ) (1), 

starting in degree 1. 

Proof. Let / be in and denote by fd- a +i x % 1 the monomials in / along 
the boundary in question. Here fd- a +i denotes a polynomial in xi and x 2 of 
degree d — a + 1. Clearly, / may be extended to an element in Q' if and only 
if we can find a polynomial gd-a{ x ii ^2) ^3 with support in Q' \ Q such that 

fd-a+l x 3 = 9d-a X 3 

(Xi + x 2 ) ■ 

The commutation relations in A imply that for any G which does not 
depend on x 3 and any C we can find f' k G A k such that 

(xt + C x 2 ) f k = f' k (an + ?i2 C x 2 ) . 

From this it follows that X\ + q a q^i 1 x 2 annihilates (Q/Q') d , as was to be 
shown. □ 

4.6 

We can summarize the discussion as follows. 



18 



Theorem 2. We have 

dcg Q 

d 3 Q^($A/(x 3 ,x 1 + q a >x 2 ) (1), (14) 

i=l 

starting in degree 1, where are the heights of the horizontal boundaries of 
Q as above. 

Proof. The Hilbert polynomial equals the constant deg Q on both sides. We 
constructed a map from the LHS to the RHS in (fl4l . Its kernel consists of 
functions / that vanish on the boundary strips along the white boundaries 
of Q and may be extended as solutions of Kf = to a horizontal strip just 
beyond the black boundaries of Q. This modified domain is a translate of 
Q(— 1) in the x% direction and our conditions on / imply / G x%Q. Thus the 
map above is an isomorphism. □ 

5 Correspondences and renormalization 
5.1 

Let Q' be obtained from Q by moving one boundary segment by one step, as 
in Sections 14.41 and 14.51 We found that the corresponding modules Q and Q' 
fit into an exact sequence of the form 

O^Q'^Q^P^O (15) 

with a point module P. Modules of the form (fT2l) form an open set in the 
moduli spaces of A-modules of rank and 

(ci,c 2 ) = (D,D(D + 3)/2) , D = degVl. 

Let M(ci,C2) denote a finite cover of this open set over which the ordering 
of the summands in (fl4"l) is chosen. 

For general C2, we define M(ci, C2) as the moduli space of A-modules with 
the same minimal resolution as the generic commutative resolution plus an 
ordering of boundary points. For example, for 

(ci, c 2 ) = (D, D(D + 3)/2 — k) , < k < D/2 , 

the corresponding modules Q" are of the form 

-> A(-2) D - k -> A k © A(-l)°- k -> Q" -> , 

generalizing (fT5|) . 
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5.2 Noncommutative shift on the Jacobian 

It is easy to see that 

dimM(ci, c 2 ) = c\ + 1 , 
and, in particular, it doesn't depend on c 2 . 

Lemma 9. The correspondence 

{(Of, Q)} C M( Cl , c 2 + 1) x M(ci, c 2 ) 
defined by (TT5T) i/ie graph of a birational map. 

Proof. It suffices to consider the commutative case, when this becomes a shift 
by P on the Jacobian of the curve Q, see Section 13.21 □ 

This means that we have an action of of a group S = Z 3ci on |_| M(ci, c 2 ) 
by birational transformations. A subgroup of the form Z 3ci_1 preserves c 2 
and acts birationally on individual components. Perhaps the title of this 
subsection is the appropriate name for this group action. 

Parallel group actions may be defined for other noncommutative surfaces. 
They turn out to encompass several previously studied discrete dynamical 
systems. This is the subject of a joint work by E. Rains and the author, the 
results of which will appear in [19]. 

5.3 

Renormalization is a central concept in mathematical physics. For any 
tileable domain Q and any m G Z >0 , the scaled domain mfl is again tileable 
and it is natural to ask how this scaling transformation affects our random 
surfaces. Equivalently, of course, one can keep fl fixed and divide the mesh 
size by m. 

To quantify the word "affects", one tries to summarize the behavior of 
random surfaces in terms of finitely many essential degrees of freedom. One 
further hopes to define an action of the group lR>o on this space extending 
the scaling transformations above. 

While this procedure is a very powerful guiding principle, in practice one 
usually has to use various approximations to make it work. This is only 
natural since one is trying to squeeze an infinite-dimensional problem into a 
finite-dimensional dynamical system. 
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5.4 



Our kind of problems are special in that they have a certain built-in finite- 
dimensionality, starting from a finite number segments that bound fl The 
quantum limit shape Q also varies in a finite-dimensional moduli space. Since 
Q encodes the essential information about the correlation functions, the 
renormalization dynamics may be considered understood once the scaling 
action on Q is determined. 

Clearly, scaling transformations are composed out of many noncommuta- 
tive shifts on the Jacobian, and more precisely we have 

Theorem 3. The scaling i— > m ■ preserves the S-orbit of Q and inter- 
twines the action of s G S with the action ofm-sES. 

5.5 

If the parameter q is adjusted simultaneously, then the limit m — ► oo becomes 
the thermodynamic limit 

mesh — > 0, logg = O(mesh) , 

which is the limit that we were planning to take all along. In this limit, 
noncommutative shifts of the Jacobian may be viewed as a perturbation of 
the commutative shift, thus joining a much-studied area of perturbations of 
integrable systems. Their analysis from this point of view will appear in [10] . 

In particular, it will be shown in [10] that, indeed the quantum limit 
shape Q is a deformation of the curve Q that defines the (classical) limit 
shape. 

6 Outlook 
6.1 

The Kasteleyn operator on Tq has an infinite-dimensional kernel and the 
Kasteleyn equation needs to be supplemented by boundary conditions in 
order to have a unique solution. By contrast, once a second difference equa- 
tion of degree d is known, the solutions form a d- dimensional linear space, 
spanned by modulated plane waves. This simple principle gives a powerful 
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way to control K -1 in the thermodynamic limit which is the basic analytic 
issue in the analysis of stepped surfaces. 

In fact, optimistically, one may expect these techniques to overcome the 
difficulties that lie in the way of proving the CLT for stepped surfaces with 
polygonal boundaries (see [7] for techniques that can handle a different sort 
of boundary conditions) as well determining the local correlations. Further, 
since polygonal boundaries are dense in the space of all boundaries, the calcu- 
lation of local correlations for them has direct implication to the classification 
of Gibbs measures (about which [20] contains a wealth of information). 



Kasteleyn theory and the formalism of (9] work for any periodic bipartite 
planar dimer. It would be very interesting to study the quantum limit shape 
in this generality. It would be also very interesting to find applications to 
difference equations other than the Kasteleyn equation, such as e.g. discrete 
Dirac equations in dimensions > 2. 



For very special kinds of boundary conditions (see e.g. [161 E] for an in- 
troduction) the g-weighted stepped surface partitions functions Z become 
generatingfunctions Zqt for the Donaldson-Thomas invariants of toric CY 
three-foldsi. 



Very generally, for any smooth projective three-fold X the expansion of 
\ogZ DT in powers of logg was conjectured to generate the Gromov-Witten 
invariants of the same three-fold X, genus by genus [TTJ. For a toric three-fold 
X, this conjecture was proven in [T2] . 

In particular, this relates the genus Gromov-Witten invariants of X 
to the leading asymptotics of logZ as logg — > 0, and hence to the limit 
shape Q. Mirror symmetry, which is certainly too complex and multifaceted 
a phenomenon to be discussed here with any precision, associates basically 
the same curve Q to X. Thus the limit shape point of view puts mirror 
symmetry on a firm probabilistic ground in this particular instance. 

The quantum limit shape Q captures the fluctuations and hence all orders 
of the expansion of log Z. This makes it a strong candidate for the as yet 

2 The customary parameter q in DT theory differs from ours by a minus sign. 



6.2 



6.3 
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mysterious higher-genus mirror of X. Moreover, the quantum limit shape 
Q, being a categorical object, might stand a better chance of generalization 
to non-toric X than the underlying box-counting. Progress in this direction 
remains both very desirable and scarce. 

6.4 

There exists a different (and, at present, conjectural) way to extract the 
higher genus GW invariants out of Q, which was proposed in [2] based on 
the diagrammatic techiques developed in the random matrix context, see in 
particular [3J. The use of noncommuting variables in a related context was 
advocated, in particular, in [5]. 

It is reasonable to expect the two approaches to converge, especially since 
various random matrix models may naturally be viewed as continuous limits 
of stepped surfaces. For example, one sees random matrices quite directly 
near the points where Q intersects the coordinate axes [TB]. There are also 
numerous parallels between random matrices and Plancherel-like measures 
on partitions, which are induced on slices of stepped surfaces [To] . 

A more ambitious goal may be to push the theory away from the Kx = 
case. As a first problem in the Kx ^ direction, one can try the equivariant 
vertex [TTJ[T2]. As a random surface model, it is very nonlocal and otherwise 
distant from what is perceived as natural in statistical mechanics. To its 
credit, it has a map, due to Nekrasov [T3], onto certain local 2-dimensinal 
lattice fermions. In contrast to Kasteleyn theory, these fermions are now 
interacting and it remains to be seen how much progress one can make in 
this more general setting. 
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